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SUMMA eee 


A method of analysis of arbitrary cascades under conditions 
of ideal fluid flow is described. The method is based on con- 
formal transformation, achieved experimentally through an elec- 
trical analogy. Tests were conducted to permit comparison of 
results obtained by this method with those of a previous ana- 
lytical solution and with NACA two dimensional cascade tests. 
Results indicate that the quick and inexpensive procedure of 
this method is applicable to engineering problems of cascade 
analysis. 

This work was conducted at the U. S. Naval Fostgraduate 


School, Monterey, California. 
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1. INTRODUCTION 


One of the first steps in the aerodynamic design of an 
axial flow turbomachine is the design of cascades which produce 
desired deflections of the fluid as it passes through the ma- 
chine. Fundamental to the solution of this problem is the de- 
termination of the two dimensional flow of an ideal fluid across 
cascades of compressor and turbine blades. An ideal fluid is 
here assumed to mean one that is homogeneous, continuous, incom- 
pressible, and inviscid. 

The development of a method which predicts the two dimen- 
Sional flow of an ideal fluid across any arbitrary cascade was 
the objective of this investigation. One méasure of cascade 
pertormance is the fluid turning angle. More basic parameters 
are the cascade coefficient K defined in Section 5.1 and the 
relative zero lift angle on defined in Section 5.1 and illus= 
trated in Fig. 17. Of primary interest was the effect of changes 
in cascade solidity and blade stagger on these indicators of cas~- 
cade performance, but a preliminary test of the ability of this 
method to predict velocity distribution around the airfoils in 
cascade is also included. 

Considerable ingenuity has been exerted in previous work 
done on this problem, and many mathematical analyses exist of 
the two dimensional potential flow past infinite cascades of 
airfoils. Some provide exact solutions for special shapes such 
as Straight blades of zero thickness. See Refs. l and 2. Some 


provide approximate solutions for particular airfoil shapes and 





cascade arrangements. See Refs. 3 and 4. Finally in regard to 
the general problem of arbitrary airfoils in cascade, there 

have been several attemyts at a mathematical solution. These 
usually involve iteration methods. See Refs. 5, 6, 7, and 8. 
The labor involved in the application of these methods is great; 
and, for some applications, even modern high speed computers 
prove inadequate. Particular difficulty is encountered with 
blades of high camber and cascade of high stagger angles. Hence 
Mesengineering application of these methods 1s Limited, 

Theoretically, the method here investigated is based upon 
the conformal transformation of an arbitrary cascade into a cir- 
cle. This is accomplished experimentally, using electrical con- 
ducting paper. The flow of an ideal fluid is analytically known 
about a circle. Through LaPlace's equation the flow can be de~ 
termined about the cascade. 

To compare the results obtained by this analog method with 
those of a previous analytical study, the airfoil in cascade 
which had been analyzed by Mutterperl and reported in Ref. 7 
was investigated. To determine the effect of solidity and stag- 
ger upon turning angle, a blade more closely approximating those 
currectly used in the gas turbine industry was desired. The 
Me 65~(15)10 airfoil section was chosen for this reason and 
because its use permitted comparison with NACA tvo dimensional 
cascade tests reported in Ref. 9. 

This work was conducted at the U. S. Naval Postgraduate 
School, Monterey, California during the 1957 - 1958 school year. 

The authors wish to express their appreciation to Professor 


M. H. Vavra and Professor T. H. Gawain of the U. S. Naval Post-= 








graduate School. This thesis project was conceived by Professor 
Vavra and is an outgrov'th of previous investigative work per- 
formed by Frofessors Vavra and Gawain. The association with 
Professor Gawain who guided and assisted in the work of this in- 


vestigation was most rewarding. 
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2. TABLE OP aawe or. 


radius of circle in the picture plane, feet 


distance from origin of sources and sinks in picture 
plane, feet 


chord length, inches or feet 
Gurtoil section coefficient of eae 


dimensionless transformation parameter, defined by 
Saud t Lon, 3..27 


SWrrehu »atiperes 
indicates imaginary part in complex number notation 


Characteristic distance used in conformal) transition — 
SLon. eet 


source strength in picture plane, feet squared per sec- 
ond 


dynamic head, pounds per square foot 

radius vector length in picture plane, feet 

tangential spacing between blades, inches or feet 

complex velocity potential function 

real axis, real plane 

imaginary axis, real plane 

complex coordinate, real plene 

potential, volts 

total current function, defined in Appendix A 

cascade coefficient, defined in Section 5.1 . Hs 
Gomeerescisctance. ohms 7 cathen deste e fron Per 

pressure coefficient = ve) 


velocity, feet per second 





to (aoe Se Oe ON OD ——~. 


angle of fluid flow, degrees 
ae absolute angle of attack 
9,6 angle of mean velocity 
X angle of zero lift 


© 


o cascade entrance angle 
6.6 


cascade exit angle 


w 


a; ideal approach angle 

cascade stagger angle, degrees 
relative angle of zero lift, degrees 
10Ge b 
complex coordinate in the picture plane 
imaginary axis, picture plane 

angular coordinate, picture plane, degrees 
specific resistance, ohms 

real axis, picture plane 

resistivity, ohm-inches 

density, slugs per cubic foot 

Solmicnay = —— 

"aspect ratio', defined by equation 3.13 
Velocity potential function 

velocity stream function 

circulation, feet squared per second 
turning angle, degrees 

vortex strength, feet squared per second 


percent potential 


Pereent current quantity 
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i upstream ot blade row 
2 downstream of blade row 
Q | axial direction 

2 ideal 

m mean 

S Stagnation 

t tangential direction 

ie leading edge 

T trailing edge 


erscripts: 





indicates maximum value 








3. ANALYSIS 


3.1 IDEAL FLUID FLOw 


If, in a two dimensional steady flov, there exists a veloc- 


ity point function such that 


V=eVvit)=V(%5y) 3.1 


the necessary and sufficient condition for this to be potential 


flov' is that 


Tx VY =0 cee 


uhich means that 
V=V@ Bue 


Pe, if further, the density of this flow be constant, continu= 


ity requires that 


Vey =o cw 


If both 3.3 and 3.4 are satisfied; then the potential is a plane 


Mermonic function such that 


Vp =0 a 


Ci 


vhich is LaPlace's equation. 
Correspondingly in this theoretical, two dimensional fluid 


mov, there is a stream function for which 


7 v =O 35 





Trajectories of constant Py are orthogonal to trajectories of 
constant wy. 

Analagous to this theoretical, tvo dimensional fluid flow 
is the flow of a steady current ina thin, flat sheet of uni- 
form conductivity. Deiails of the analogy are contained in Ap- 
pendix A. On the sheet the electric potential E, analagous to 


QP , can be represented by 
vie re. 3.7 
Sematarly, introduction of a total current function 1 results an 
wile 3.8 


Trajectories of constant E are orthogonal to trajectories of 
constant 1. 


3.2 CONFGRMAL TRANSFORMATION 


Various flows of an ideal fluid are analytically known a- 
bout certain simple domains such as a circle. LaPlace's equa- 
tion, which characterizes such a perfect fluid flow, remains in- 
Wariant for any conformal transformation. Hence, if an arbitrary 
shape can be conformally transformed into a suitable circle, the 
flow of an ideal fluid about the arbitrary shape can be deter- 
mined. 

The theory of conformal transformation with the aid of an 
electric tank is discussed in Ref mer In general the whole re- 
gion between two arbitrarily specified closedcurves is not trans- 
formable into the whole region between two other arbitrarily 


Specified closed curves. Transformation can be accomplished, 





however, if there exists a function, in the picture plane, of 
the real plane complex variable 2 such that for every assigned 
value of the variable the function has a definite value or sys- 
tem of values and also a definite derivative. 

The arbitrary shape transformed in this investigation con- 
Sists of an infinite cascade of equally spaced compressor 
blades of the same shape and orientation. 

Imagine a conducting sheet of uniform conductivity, infin- 
ite in extent in two dimensions, Fig. 1, in which the blades 
Sire areas of infinite resistance (cutouts from the sheet). With 
a potential E* imposed on a line parallel to the axis of the 
cascade and located infinitely far ahead of the cascade and a 
zero potential imposed on a similar line located infinitely far 
behind the cascade, corresponding points on each of the blades 
have equal potential. Immediately in front of the cascade a 
line connecting points of equal potential exhibits a shape that 
is strongly influenced by the geometry of the cascade; but with 
increasing separation from the cascade the influence diminishes 
rapidly, and at finite distances ahead of and behind the cascade 
these lines of constant potential are virtually straight and 
parallel to the axis of the cascade. Hence, the conditions at 
infinity may be closely approximated by imposing a voltage on 
the sheet through straight bus bars parallel to the axis of the 
cascade and a finite distance away. 

No generality is lost if a single strip of width s which in- 
cludes the space between two blade cutouts is considered, provid- 
ed that the boundaries or edges of the strip are lines of con- 


Stant 1. Let the boundaries be chosen to coincide with portions 





cr 


of the blade profiles as shown in Fig. Z. 


For the strip of conducting sheet 
R¥iL*=E* S.3 
For points on this strip (from Appendix A) 
BE +épl = $(2) 3.10 
Mividing both sides of 3.10 by 3.9 


fil2l= Esti ba sEtid 3.11 


Sill be constant along the edge of the strip and E will vary 
zrom one on the front bus bar to zero on the after bus bar. Any 
point on the strip can be located by coordinates x, y; but © and 
wk are continuous functions on the strip and are definite at any 
point. Optionally, then, any point may be designated by its co- 
ordinates or by the values of & and aN which characterize it. 
In symbols: 
€ =E(x,y) x =x(C, H) 


and conversely Sake 


ded(x,y) 7 =y(6,&) 


On the front part of the profile cut-out © will reach a maxi-~ 
mum value at a point 1, and on the trailing part a minimum 
Marve at a point 2. Fig. 3 shows a typical plot of percent po- 
tential S vs. percent chord. The quantity € -€, 1S Gne Wen 
the parameters used to accomplish transformation. 

In the region near the bus bars where lines of constant & 


are straight and parallel to the cascade, imagine two lines of 


10 





constant & <a and or such thaté,—c,=C, - 6, . see Figwmae 


Imagine that the width of the conducting strip, s, is divided 
mito 1), equal parts, each measuring 3 This same increment 
| 


of length is used to measure along the distance Bl between Og 
_——- case | _ og 
and om dividing this distance into n= Tn 


Sonstant a drawn at these stations are separated by equal digs 





parts. Lines of 


tances and represent equal increments AE, . For convenience it 
Memassumed in Fig. 4 that nN, and i, are integers. Lines of 
constant é and constant ok drawn at these measured stations 
Meem an Orthogonal grid of NM, Squares over this portion oF 
eae strip. 

A similar conceptual procedure is valid at any position in 
the strip. Lines of constant & may be drawn such that 
é, —é, S E, —C, . Let the distance between them, as meactimed 
along a line of constant 1, be divided into 1, segments of equal 
potential drop, 46 » The distance Bere icici the edges of the 
Strip, as measured along a line of constant & ,» can then be die 
Vided into 7, segments. A grid of orthogonal curvilinear co- 
ordinates results and the region consists of the same number 
NN, of curvilinear squares as before. 


Let there be a eee tT such that 


(earn 


ae sce, Sere 


In the regions of the bus bars a varies linearly with 1; so 


(i) = se 3.14 


t then reduces to the ratio & which might be considered an 


i 


"aspect ratio! for the gridded section. lt represents the ra- 


ia 





tio of the number of squares along the Strip to the number of 
Squares across the strip. The quantity represents a unique and 
invariant characteristic of the geometry, and must be maintained 
in mapping to achieve conformal transformation. 

Consider the denominator of the quantity @. From the anal- 


ogy developed in Appendix A, at any point, 


| SE 
Vix = 2 3.15 
‘7 x 
In the region near the bus bars where ay is constant across the 
width of the strip 
Uae Pa 3.16 


ye 


Multiplying both sides by s, the width of the strip, 


7. oe (E.- Ee) (+ . cf 


Ds is the total current flowing in the strip between bus bars. 
The right side of the equation might be labeled the specific 
voltage drop of the strip. It represents the voltage drop be- 
tween two equipotential lines which form a Square with lines of 
constant I that are s distance apart. Equations 3.15, 3.16, and 
3.17 may be used in non-dimensional form if each is divided by 
equation 3.9. 

Now consider the 8 Or picture plane of Fig. 4 with sources 
of equal strength m on the negative, real axis at -b and -= 
and sinks of the same strength on the positive, real axis at +b 
and + . These sources and sinks create a complex function 


in the 5 plane. See Appendix B: 


ee 





[ 
eee abr S ~ Dr 3? 
ur =o : b-® + bn bt zy, te ee 3.18 


Any point in the picture plane can be located by coordinates Ns 
& or by the values of PP and Y which characterize the point. 
symbolically, 


P= P(EsX) E=E(P; y) 


and conversely 3.8 


ye wlEs2) | n= n(P, y) 


The circle [5 Sa 1S a line of constane Y ; MOuUGhL. 1S etme 





conformal transformation of the strip of the real plane external 
to the blade cut-out onto the area of the picture plane outside 
the circle so that the cascade blade outline becomes the circum- 
ference of the circle. p is a minimum at a point 1 and a maxi- 


mum at a point 2. The difference in potential L-P is 


_m coshy +] 
R-G= Fn Se aa 


Lae coh*F¥+/ _~ Z-P aaa 


a 


1T coshy—-| mM 
See Appendix B. It will be shown that the right side of eq. 3.21 
is the corresponding form of @ for the picture plane. 

With reference to the strength term m in the denominator of 
eq. 3.21, recall that in two dimensional, unidirectional flow of 
an ideal fluid the uniform quantity rate of flow between stream- 
lines h distance apart is Vh, Fig. 5. Correspondingly for a 
point source from which flow emanates, a uniform quantity rate 
of flow mroeces any Simply connected curve enclosing the source. 


The velocity at any point is the gradient of p ; “chia bees Vac 


ie 





The change of p between two DP countours which, with two 
Streamlines h distance apart, form a curvilinear square is 
<Eh ». ~his quantity might be dubbed the specific P drooe 
analogous to specific voltage drop, but for fluid flow consid- 
erations it is designated a quantity rate of flow, here repre- 
sented by m. 

To achieve conformal transformation between the z and 5 


planes, the following condition is imposed 
P=k,- k2€ 3.22 
where K, and K, are arbitrary constants. 


It follows that Q-P =k,(6,-E,). 3.23 


Further, let M=K25 (2c) | ° 3524 


Then = KR fe) 33.29 
28 aes | 
3X Joe 

The consequence of equal @ in the real and picture planes is 
this: for any area in the real plane bounded by lines of con- 
Stant fe, and a there is an area in the picture plane bounded by 
corresponding lines of constant Pp and WwW » The 'aspect ratios' 
of the areas are the same. See Fig. 4. 


For every point Ss ; wk in the real plane there is a cor- 


Ao 
Az 





responding point DP, Ww in the picture plane. can be meas- 


ured in the real plane, and a corresponding a can be calcu- 
lated in the picture plane. These conditions are sufficient to 
establish the existence of a function of a complex variable which 


relates the two planes conformally. An analytical expression of 


this function is not determined nor required. 


14 


Definition of @ in the 5 plane uniquely determines 
the location of the sources and sinks. Eq. 3.21 is solved giv- 
ing 

We. 
coshy = © +t 3526 
err | 

Further, the periodicity of the mapping function is such that 
each strip of width s in the infinite cascade is mapped onto 
this same area in the picture plane; so the result is truly a 
map of the infinite cascade. Fig. 6 indicates the location of 
corresponding points on the two planes. 

Reference to eqs. 3.23, 3.24 and 3.25 indicates that only 
three quantities need be measured in the real plane to accom- 
plish the transformation. These are 6-6, 60). S 
Other measurable quantities might perhaps be used, but the a- 
bove listed parameters seem to be the ones most naturally, con- 
veniently, and accurately obtainable. They are sufficient to 
achieve transformation. 

The position of a point on the blade profile such as P in 
Fig. 6 may be expressed in terms of the measured value ofae 
which characterizes it. Similarly the position of a point on 
the circle such as P' in Fig. 6 may be expressed in terms of the 
calculated value of YP which characterizes it. Points P onvege 
blade and P' on the circle which correspond have identical val- 


ues of the parameter f, where f is defined by tie expression 


otal P-L 
€-67 ee 5-9 3.27 


Each such point P' on the circle also has an angular co- 


ordinate © as shown in Fig. 7. The coordinate OG is related to 


les 





the parameter f by the expression 


J | 
= cos lcoshy 
Gr cos Sno tt, ) 3.98 


This formula is derived in Appendix B. 


3.3 THE FLOW BER ie. 


Once the transformation has been completed the flow in the 
picture plane can be modified to correspond to various desired 
flow conditions in the real plane. This involves variation in 
the direction of the mean flow and variation in the circulation 
around the blade to satisfy the Kutta condition. 

The direction of the mean flow may be changed by super- 
position of a set of vortices /J\ at the Singularities in the 
picture plane as shown in Fig. 7. The circulation around the 
airfoil can be modified by superposition of a second set of 
vortices aie as indicated in the same figure. 

The Kutta condition states that the rear Stagnation point 
of the fluid flow about an airfoil with a sharp trailing edge 
occurs at the sharp trailing. edge. Therefore the rear stagna- 
tion point in the picture plane must occur at the point on the 
circle which corresponds to the trailing edge. This determines 
the strength of the circulation that must be imposed. 

At the desired angle of attack and with the Kutta condition 
Satisfied, the velocity Vy at every point on the circle may be 
calculated. Since the conformal transformation has established 
the correspondence between points on the blade in the real plane 


and on the circle in the picture plane, the velocity V2 at each 


ic 





point on the blade is obtainable by means of the relation 


Vj 
had El 3.99 

, en 
35 ap (A=£-) 3.80 











It may further be shown that d2 = 


ae 
aA 7 
It should be noted that when transformation has been accomplished 


where is measured along the blade profile. See Appendix B. 


velocities of corresponding points on the blade and circle are 


dz 


related by the factor , which is a constant at any given 








point and is not a function of angle of attack or circulation im- 
posed upon the circle in the picture plane. 

The detailed analytical development of the ideal fluid flow 
in the picture plane is contained in Appendix C. What follows 
here is a descriptive summary of the procedure. The sources m 
placed in the picture plane to accomplish the conformal trans- 
formation create a ‘straight through! fluid flow. This flow 
causes the front and rear stagnation points to occur on the cir- 
cle at the ends of a diameter which lies along the real axis. 
meewpoints 1 and 2 in Fig. 7. 

To develop an angle of attack of the fluid with the zero 
direction of the real axis, vortices of atveaerne are imposed 
upon the four point sources. For positive approach angles, coun- 
ter-clockwise vortices of equal strength are placed at the source 
Mee —6 and at the sink at&= = . Clockwise vortices of 


the same strength are placed at the other source and sink. The 


mean approach angle then becomes: 


ot Ve 


Arm = tan == 3. en 


ae 





The front and rear stagnation points still occur at the ends of 
a diameter of the circle but the diameter makes an angle with 


the real axis 


Os = tant’ (tan Xp tanh ®t) | CHE 


The flow is undeflected as it passes around the circle; so no 
net circulation or lift is present. 

To produce deflection and lift, vortices each of strength 
ai and of the same direction are placed at the source and sink 
inside the circle. The strength and the direction of these 
vortices are those required to satisfy the Kutta condition; 
i.e. the rear stagnation point must be brought to the location 
on the circle which corresponds through mapping to the trailing 
edge of the airfoil. Vortices of the opposite direction must 
be placed at the source and sink outside the circle to insure 
that the stagnation streamline remains a circle of radius a. 

With sources and sinks m, and with vortices me and £ 
in place, the picture plane is complete. The velocity at any 


point on the circle is 


Va = 2m cosh ¥sin@ }+ sinh ¥cos I” { sinhev 
Vv oe O7~ 7a \cosh2Zt-sinch} pa \coshay—seneB) 27a\cosh27-sMzO ] . 2. 


feo cetermined by m= = s(Va) = s(- x) 3.34 


IN 


aa is determined by eq. 3.31 once Ky, 18 arbitrarily chosen. 
To satisfy the Kutta condition, let 6=6,, and Vo =O. The val- 
ue of QO, is determined by substituting the value of $7 in 


Spe. 25. SOlving eq. 3.33 then gives 
i cosGr sin Gr 
eer coshe — sinh? 3.35 
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The picture plane is thus completely determined. See Fig. 7. 
For subsequent work, formulas for particular quantities 
relating to the flow direction are developed. The cascade en- 


trance angle is 
A= fan [tank + gi oe ee 


The cascade leaving angle is 


The turning angle of the fluid as it passes through the cascade 


is 
GQ =XK,-& 3.38 


Of aerodynamic interest is the ideal angle of attack. When the 
front stagnation point occurs at the point of the circle corres- 


ponding to the leading edge of the airfoil 


= = } rm sinGr — sin, 
Ba= i tanh sorte ea oie 


The lift on any blade is developed from the elementary relation 
| = PVmT" =X, 141 3.40 


iS is the density of the fluid and A is the projected area of 
the section along the chord. The subscript X appearing in the 
above listed equations indicates that the coefficient of lift 
and the dynamic head qd may be based upon the cascade entrance 
velocity V , the cascade leaving velocity Vp , or the cascade 


mean velocity ve Each of the three coefficients of lift has 


its reasons for being used in the analysis of cascade performance. 
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fy , based upon Ve » permits an aerodynamic tie-in with single 
airfoil theory, <g » based upon V; is useful in cascade design 
because conditions upstream of the cascade are generally known. 
xf, based upon We , makes allovance for the deceleration of 
the flow as it passes through a compressor cascade. Ref. ll 
states that acutal cascade tests have shown that the drag-lift 


ratio reaches a minimum at an approximately constant value of 


Aa = O.58 ® 
In terms of the quantities determined by the previous e- 


@Meur0ns, the coefficient of lift is 


_ 2 608 Kx 
Abe = SE (tan, — tanog) 3.41 
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4. EQUIPMENT 


The equipment consisted of the specially designed apparatus 
used to hold the two dimensional, conducting paper, with appurte- 
nances to impose the potential and complete the circuit. Associ- 
ated with this piece of equipment were various electrical in- 
struments necessary for potential supply and meaSurements. A 
list of the instruments with specifications, where applicable, 

5 Given in Table 1. 

A complete test setup is shown in Fig. 8. The figure shows 
the conducting paper, the working surface and the associated e- 
lectrical instruments. The remaining item is the profile tem- 


plate used as a guide in cutting the conducting paper. 


4.1 THE WORKING SURFACE 


Fig. 9 shows the original working surface which existed at 
the start of the investigation. 1t was constructed of 3/4 inch 
plywood overlaid with 3/16 inch Masonite (rough surface toward 
the plywood) and was fitted with Birmingham Gage Number 25 copper 
bus bars, one inch wide and five feet long. These bus bars were 
positioned on the pile parallel and 29.5 inches apart at their 
inside edges. They served as the electrodes for imposing a po- 
mmetal difference across the paper cascade. 

The paper was placed over the bus bars with its conducting 
(graphite) surface toward the Masonite. It was clamped as 
Bown in Fig. 10. This method of clamping provided a tight, 


line contact between the bus bars and the paper. 
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This setup left the ends of the paper as unknown boundaries 
which vere streamlines of the flov. It was necessary to uSe an 
iteration process to determine these boundaries. A first approx- 
imation was made by eye and the paper cut accordingly. The suc- 
ceeding approximations were then obtained by flux plotting a 
streamline from experimentally determined constant potential 
lines. These plots were made near the middle of the cascade to 
minimize effects of the previous boundaries. The cascade ends 
were then cut to this nevly determined streamline. This process 
of refinement was continued until the correction required became 
as small as the errors of measurement. 

In order to by-pass the time and labor consumed in obtaining 
these boundaries a second table was constructed. This table 
Shown in Fig. ll, was similar to the original in its method of 
connecting electrodes to the conducting paper. The contact be- 
tween the electrodes and the paper was obtained by the same man- 
ner of clamping. The copper bus bars were widened to 32 inches 
however, and maintained parallel to + 0.003 inches in a three 
foot length. 

The table vias lengthened to six feet overall with a foot 
and a half of each end hinged to the main, three foot section. 
“With this provision, the paper was clamped to the table as in 
Figs. lO and 12 and the ends were folded down under the main 
section thereby bringing the ends of the paper Poeace together 
in a physical and therefore electrical contact. The cascade 
then had no beginning or end and was, in effect, infinite. 

Since cascades of varying solidity required various lengths 


of paper in order to come to contact at the joint, it vas neces~ 
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sary only to lengthen or shorten the bight of the paper at the 
hinges when securing the paper to the table prior to folding. 
The method of obtaining the electrical connection between 
the ends of the paper cascade is shown in Fig. 13. This method 
has the disadvantage that the operator cannot see the actin 
connection and must rely upon the trueness of the hinges, ma- 
chined parts, etc. To provide sufficient contact pressure a- 


long the joint, C-clamps were used outside the aluminum bars. 


4.2 THE CONDUCTING MEDIUM 


Type "L"' Teledeltos paper was used as the conducting medi- 
um throughout the investigation. The specific resistance of 
this paper is given by the manufacturer as 4000 ohms, with an 
accuracy of + 3 to + 8 percent with the lower value of resist- 
ance along the length of the roll and the higher value across 
the width of the roll. The quantity specific resistance, is 
defined and explained in Appendix D. The corresponding specif- 
ic resistance of the copper bus bars is clearly negligible in 
comparison with that of the paper. At the hinged ends of the 
table, and near the electrical joint on the under surface, where 
the bus bars are necessarily interrupted, the required conduc- 
tivity was secured by painting the paper on its conducting side 
with a special silver paint supplied by the manufacturer. The 
Specific resistance of such painted areas iS given by the manu- 
facturer as about one to four ohms. While this value is consid- 
erably higher than that for the copper bus bars, it is still 
negligible in comparison with that of the paper. Thus the con- 


ductivity of either the bus bars or of the painted areas is es- 
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sentially infinite in comparison with that of the paper. Each 
test setup required a new piece of paper making it a consumable 


item in the investigation. 


4.3 THE BLADE TEMPLATES 


The following profiles were selected for investigation: 

a. The Mutterper] profile, ordinates presented in 
Ref. 7 and listed in Table 1l 

b. A typical compressor profile, the NACA 65-(15)10 
section, ordinates presented in Ref. 9 and listed 
in Table III 

Coordinates ot these profiles, scaled to the chord sizes used 
in this investigation, are presented in Tables IV and V. 

In order to transpose the coordinates to the paper in the 
proper geometry and to maintain accuracy, templates were used as 
cutting guides. They were made of 1/16 inch, oil hardening flat 
Stock which was left in the spheroidized condition instead of 
being hardened, thus avoiding distortions. 

Tolerance from the tabulated coordinates was held to +.005 
inches. Each template was checked dimensionally at fifteen dia- 
meters magnification on the Bausch and Lomb Contour Measuring 
Projector, Model 33-12-01, at the Supersonic Wind Tunnel of the 
U. S. Naval Postgraduate School. 

To position the templates for cutting on the cascade layout 
lines, three or more locator holes were drilled in each. These 
holes were spaced along a line parallel to the chord and offset 
from it. One of the holes was located at midchord. These pio- 


files are shown in Fig. 14. 
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4.4 ELECTRICAL CIRCUITRY AND INSTRUMENTS 


The electrical circuitof the continuous cascade table and of 
the associated instruments for measuring percentage potential 
is shown in Fig. 15. The power supply used was a Heathkit Bat- 
tery Eliminator, Model BE-4. The voltage output was maintained 
at about ll volts in order to give a quicker and more definite 
response on the null meter. The null meter used was a General 
Electric DC Micro Ammeter, Type DO-71, Model A292AA1 of 10 mi- 
Broamperes full scale deflection. At a 100 percent potential 
equal to about 1l volts, ©.01 percent potential was the minimum 
discernable change on this meter. The voltage divider, manu- 
factured by Electro-Measurements Inc., was a hexadecdde type; 
Model RV622. Its linearity as cited by the manufacturer was 
+ 0.0025 percent with an input resistance of 10,000 ohms. As 
mentioned above, 0.01 percent potential was the minimum detect- 
able change and therefore only the four higher decades of the 
voltage divider were used. The test probe or stylusS was Sup- 
plied with the General Electric Analog Field Plotter Kit. It 
was fitted with a shallow tapered needle from a drafting in- 
Strument. 

The approximate derivative of the mapping function could 
be determined by graphical means or by direct measurement. The 
latter method was used. To accomplish this determination, the 
potential difference over a small distance was measured by a 
double electrical probe. This probe was adapted from a small 


bow compass by placing a second, insulated needle in the lead 





receptacle and soldering both needles to hookup wires. The com= 
pass was adjusted to the smallest needle gap possible and locked 
in position. This gap measured 0.04048 inches on the measuring 
projector previously mentioned. The entire compass was finally 
covered with electrical insulating tape. 

Fig. 16 shows the circuit schematic for determining the ap- 
proximate derivative of the mapping function. The hookup wires 
from the probe were lead to a Leeds and Northrup Portable Pre- 
cision Potentiometer No. 8662, where they were connected to the 
E.M.F. binding posts. This instrument has a range of from O to 
65 millivolts. The power supply was two six volt batteries con- 
nected in series to give approximately 12 volts. Fractically all 
Stations on the profile gave voltage increments across the double 
probe which were within the range of the potentiometer at a bat- 
tery supply of about 12 volts. The voltage across the cascade 
(100 percent potential) was measured to the nearest 0.01 volte 


On a Weston DC Voltmeter, Model 45, having 0-15, 0-150 volts full 


scale deflection. 
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5. TEST PROGRAM AND PROCEDURE 


The test program was designed to show the effects of various 
geometric parameters upon the performance of a cascade. The param- 
eters which define the geometry of two dimensional, rectilinear, 
infinite, cascades are the blade profile, the stagger angle p, 
and the solidity c/s. The quantities 6 , ¢ and s are illustrated 


meerag. 17. 


5.1 DEFINITION OF TERMS 


The direction of the incoming flow may be specified in terms 
of the approach angle & , which is measured from a line perpen- 
Mmeutar to the cascade. 

An alternative parameter which is sometimes useful is the 
absolute angle of attack &q, which is defined as the angle be- 
tween the mean flow direction and the direction of zero lift, 
Axis 1. The mean flow angle Oy, , depends both upon the inflow 
and outflow angles oj and &, . The angle of zero lift oy, ; 
of the cascade is not initially known, since it depends on the 
entire cascade geometry in some complex manner, but it may be 
found from the results of the test measurements in any given 
case. Once &, has been determined, a convenient parameter for 
further analysis of the cascade is the relative angle of zero 
fet t Ei. » which is defined as the angle between the chord and 
Axis 1, the direction of zero lift of the profile in cascade. 

These various relations are shown in Fig. l7. 


The performance of a cascade may be expressed in terms of 
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any convenient parameter which relates the inflow and outflow 
directions. Among these are the turning angle @ and the lift 
coefficients ,; Kp, or C4, . The turning angle @ , repre- 
sents the net change in direction of flow from inlet to exit. 
The lift coefficients Cy , &g Andee, represent thew. 
tio of the mean effective pressure difference across the blade 
to some reference dynamic pressure. They differ only in the ref- 
erence condition used: Cp uses the mean velocity head, while 
Kg, and «< g. employ the iniet and outlet velocity heads re- 
Spectively. All three of these coefficients are closei, relat- 
ed to the turning angle, and the choice of parameters is large- 
ly a matter of convenience, although for a real fluid the quan- 
ity Kp, appears to be an approximate criterion of blade 
stall. 


It can be shown theoretically that for an isolated airfoil 


cp = 2 Sind ayy, 1 
Similarily, for a cascade of airfoils, it can be shown that 


pop = 27K Sind a 


where &q = Am—ol,, and K is the cascade coefficient which is a 
constant determined by cascade geometry. 

There are two basic parameters, then, which completely de- 
termine the performance of a given cascade in ideal fluid flow. 
These are K and P,,. . They have been calculated for the cas- 


cades tested and are reported. 
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5.2) THE TEST Peer 


The method used in this investigation requires that a sep- 
arate set of measurements be made for each cascade configuration 
tested, whereupon the effect of variations in inflow direction 
can be determined analytically. This greatly reduces the amount 
of testing that is required; nevertheless, the testing procedure 
was necessarily painstaking, and due to the developmental nature 
of this investigation only a limited number of tests could be 
undertaken and analyzed in the time available for this study. 

The various combinations of blade profile, solidity and 
Stagger angle which were investigated are summarized in Table V1 
and discussed below. 

The airfoil shape presented by Mutterperl in Ref. 7 was se- 
lected for the first series of tests. Published data on this 
profile had been determined analytically. Tests on the profile 
in this investigation included solidities of 1.0 and 1.5, both 
at a stagger angle of 45%. The latter case was not covered by 
Mutterperl in his report. 

The second series of tests was designed around a typical 
compressor profile at combinations of stagger angle and solidity 
Characteristic of practical compressor geometry (based on the 
authors! experience with the Gas Turbine Division of General 
Electric during the 1957 Summer Field Trip from the U. S. Naval 
Postgraduate School). 

Two groups of tests were conducted using an NACA 65-(15)10 
profile at stagger angles of 20° 47 20" 30°, and 35°. One group 


used the six inch chord at a solidity © = 1.0 while the other 





used the twelve inch chord at 0 = 1.5. At least one check 
point for each test was available from the experimental data 
presented in Ref. 9. 

One of the objectives of this investigation was to demon- 
strate an extension of the method namely, the determination of 
the velocity distribution around a blade in cascade. Data for 
this demonstration was gathered during the test of NACA 65-(15) 
10 (12" chord) at © = 1.5 and ® = 35°. Selectiono this test 
provided the Opportunity to compare results with those publish- 


ed in Ref. 9. 


5. oS PROCEDURE 


The procedure is perhaps best illustrated by following a 
typical cascade setup through to data collection. 

The first step of the procedure was preparation of the con- 
ducting paper for the table. The paper was supplied in rolls, 
34 inches wide and about 20 feet long. The working surface had 
been designed with the electrodes positioned to accommodate 
this width so it was necessary only to trim the paper to length. 
The length required varied with solidity for a given blade chord. 
For @ = 1.0 and a six inch chord, 88 inches of paper were used. 
For the tests with 0 = 1.5 and the twelve inch chord, 92 in- 
ches were required. 

The paper was laid out flat on a smooth table and fastened 
without stretching by masking tape at intervals along the edges. 
The black (non-conducting) surface was exposed as the drawing 


Surface for layout lines. 
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Fig. 18 shows a sketch of the layout lines on a typical 
piece of paper. Fig. 19 shows the template outline on one set 
of these lines. Included in this latter sketch are the refer- 
ence lines near the leading and trailing edges which were used 
to locate data points. 

The cascade center line was drawn first. The best results 
in keeping this line straight occured when using a chalkline of 
sewing thread and blackboard chalk. The line so obtained was 
then penciled in using a six foot straight edge. Exact center- 
ing of this line on the paper was not necessary as the copper 
electrodes, being 32 inches apart, allowed some deviation. 

The cascade center line was used as the reference for lay- 
out of the selected spacing and stagger angles. Care was taken 
to locate the profiles so that the electrical joint would occur 
midway between adjacent blades along this reference line. In- 
Stances of large stagger angle, where the seam passed through 
corresponding sections of the cutouts on each end of the paper, 
could be handled by the joining device. 

A drafting machine was utilized to construct the angular 
stagger and midchord lines shown in Figs. 18 and 19. 

mi order to maintain accuracy in locating data points, ref=> 
erence lines were drawn on the three center profiles of the pa- 
per cascade. These lines were construgted perpendicular to the 
Chord and extended to intersect the contour of the profile cut- 
outs as shown in Fig. 19. The interval along the chord of these 
lines was 1/40 inch for the six inch chord and 1/20 inch for the 
twelve inch chord. A jewelers' eyepiece, and 6H pencil, and the 


scale of the drafting machine arm were used to construct these 
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lines. Only part of the profile was lined in this manner. Fif- 
teen lines were drawn at the trailing edge in all tests. The num- 
ber and position of the lines near the leading edge varied with 
mee location of the front stagnation point. 

The paper was then removed to a Masonite surface for cutting 
out the profiles. In this operation the template was located by 
means of the holes provided. Cutting proceeded from the tail 
and from the nose toward the midchord where the cuts met. The 
knife used for cutting was an "Xacto' using the blade having the 
Straight cutting edge. The blade was onedunte a shape illus- 
trated in Fig. 20 and stropped on a Masonite surface to remove 
all Brees . In cutting, it was imperative to hold the knife so 
the cutting edge was at an angle of about 45° to the paper and 
to guide it in toward the template. All other methods of hold- 
ing and cutting either did not adhere to the profile shape or 
produced jagged edges. | 

The ends of the paper were then folded to fit over the 
corner of the aluminum bars at the joint of the table. The fold 
line was marked on the non-conducting (black) surface so the pa- 
per was first folded to keep this line exposed as a guide. After 
this fold had been sharply creased, it was reversed to expose the 
conducting surface. | 

Painting the bus bar extensions was the next step in the 
preparation of the paper cascade. To provide guidelines for this 
painting, the paper was marked while in its proper position on 
the working surface. A center line had been scribed on the table 
during its construction. Bringing the cascade center line into 


cOincidence with the scribed line properly positioned the paper 
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in one direction. Positioning in the other direction was ac- 
complished by fitting the folds at each end of the paper to the 
right angle corner on their respective parts of the joining de- 
vice. This is shown in cross section in Fig. 13, with the paper 
joined. The end sections of the paper were held by suitable 
weights on each end section of the table while the middle of the 
paper was moved along the scribe line of the main table section 
until the bights of the paper over the hinges were equalized. 
Fig. 12 shows the paper clamped in this position. The guides 
for the bus bar extensions were then marked on the paper's con- 
ducting surface by pressing it down over the sharp inner edge of 
the ends of the copper bus bars. 

The paper was then removed to the drawing surface and layed 
out flat with its conducting surface up. At each end of tlh pa- 
per the painting guides were pencil lines which had been extend- 
ed along the bus bar impressions to the ends of the paper. 
Painting guides for the bights of the paper were obtained by pen- 
cil lines connecting the bus bar impressions at the ends of each 
bight. Scotch Brand Electric Tape was used as a masking tape. 
It was placed along the above lines on the side away from the 
paper edge and extended beyond the bus bar marks. The edges of 
the paper were then painted with the silver conducting paint a- 
long each of the masked areas. When the paint had set but not 
dried, the masking tape was removed. Unless this tape had been 
dusted lightly with blackboard chalk or some other dust, the 
graphite conducting medium was liable to pull away from the pa- 
per during tape removal. 


At this point various reinforcing sections of tape were ap- 
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plied to the paper. This was especially needed when using the 
twelve inch profile cutouts. 

After the paint had dried, the paper was aligned on the ex- 
tended table, non-conducting surface up as when marking, and 
held. The electrical tape was used in small pieces to hold the 
paper in position by placing it ion the paper edges behind the 
contact line of the electrodes and over blade cutouts not used 
for data. The clamps were then applied as shown in Figs. 10 and 
m2. 

The ends of the paper sticking through the slots of the a- 
luminum bars were cut off to leave about 3/8 inch of paper ex- 
tending beyond the fold. A 3/8 inch wide strip of electrical 
tape applied to the cut edge of these ends was used to hold the 
fold over the corner of its respective bar. The relief in the 
movable bar gave room for the tape and at the joint itself only 
the two layers of paper, conducting sides together, remained. 

The table was turned over and the ends folded toward each 
other. The folds of the paper were joined simultaneously and a 
Support under the two ends prevented their moving through their 
closest point of approach. The ends of the table were secured 
against the middle by bolts. Insulating paper was then slipped 
between the conducting paper and aluminum bar and pushed to the 
joint to prevent the curled cutout edges from contacting the a- 
luminum bars. C-clamps were applied to press the bars together 
and make a tight contact between the folds on the ends of the 
paper. 

The table was turned over and the circuit connected as in 


meg. 15. 
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The data points were numbered from leading edge to trailing 
edge. Stations every 1/& inch along the chord on the six inch 
profile were numbered consecutively. Data points between these 
lines were lettered a, b, c, and d. The data consisted of the 
value of © at the leading edge, the points near a , and at 
the fifteen upper and lower points near the trailing edge. To 
obtain a. , a line perpendicular to the electrodes and the 
cascade center line was drawn from the electrode ahead of the 
cascade. It was marked at intervals of one inch for five or six 
inches from the bus bar. The value of oS was recorded for each 
of these points. All values of © were read and recorded to 


0.01 percent potential. 


aG 
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When obtaining the derivative approximation it was 
necessary to use the circuit of Fig. 16. Percent potential 
could be obtained simultaneously by connecting the voltage di- 
vider with null meter and probe across the battery supply. 

The data in this test consisted of G at various points a- 
round the profile edge, corresponding values of AE in millivolts 
Obtained from the precision potentiometer and total potential E* 
in volts. The double probe was positioned so that its interval 
between probes was centered about the point where © was obtain- 
ed and with its needles just making contact at the edge of the 


profile. The value of De was read and recorded to the nearest 


foie millivolt. 





6. DEVELOPMENT opera wei. 


A sample of the data taken during these tests for the deter- 
mination of turning angles is presented in Table VII. A sample 
of the data taken for the determination of velocity is presented 
in Table VIII. Table IX contains the input quantities used for 
the determination of turning angles herein reported. 

Results of tests on the cascade analyzed in Ref. 7 are pre- 
sented in Fig. 2l as a curve of coefficient of lift versus mean 
velocity direction. 

Results of tests made on the NACA 65-(15)10 airfoil are pre- 
sented in Figs. 22 through 3l. In Figs. 22 and 25 the values of 
turning angle for various blade stagger angles are plotted versus 
cascade entrance angle with lines of constant coefficient of lift 
ep meossplotted. Fig. 22 is for a solidity of 1.0 and FigvgZze 
MemtOr a solidity of 1.5. Figs. 23 and 26 contain, for a solid- 
ity of 1.0 and 1.5 respectively, curves of coefficient of lift oy 
for various blade stagger angles, as a function of cascade en- 
trance angle. Figs. 24 and 27 show, for a solidity of 1.0 and 1.5 
respectively, curves of coefficient of lift fp, for various blade 
Stagger angles, as a function of cascade entrance angle. Fig. 28 
is a plot of cascade coefficient K versus blade stagger angle © 
Fig. 29 is a plot of relative zero lift angle OL. versus blade 
stagger angle 6 

Figs. 30 and 31 contain the results of the test of this meth- 
od to predict velocity distribution about the NACA 65-(15)10 air- 
foil in cascade at stagger angles of 35 and solidity of 1.5. The 


= \ 
results are plotted as (2) versus percentage distance along the 
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chord. 

The recorded data necessary to determine turning angles 
and coefficients of lift for a given cascade (blade, stagger, 
and solidity) consisted in part of readings of the percentage 
voltage at several stations of known location in the vicinity 
of je, 4 oe , and Os ». oee Table VII. Data were recorded for 
at least three blades in each caScade test, but final values of 
this report are based on the median blade for each case. These 
percentage readings near a, and ‘os could be silatiees versus 
distance along the chord to permit utilization of several volt- 
age readings to insure greater accuracy in the determination of 
the values of oO, Aner . See Fig. 32. 

Bonding om was a more difficult problem. Fig. 3 indicates 
that the curve BE © versus percentage of the chord becomes tan- 
gent to the line representing 100% chord at the trailing edge. 
Determination of ome the point of tangency, by experimental 
measurement meant reliance upon a Single potential reading be- 
cause of the physical nature of the very sharp trailing edge. 

The method used in this investigation to determine ee more 
accurately involved obtaining an analytical curve of eo near the 
trailing edge of the airfoil. The form of the curve was that of 
potential flow around a straight sided wedge. The curve was ob- 
tained by fitting this form to several experimental points in a 
finite but small distance near the trailing edge. Curve fitting 
was accomplished by means of the method of least squares. The 
point of tangency or was then determined analytically. Appen- 


dix E contains the details of this curve fitting procedure. 
Xe 
ox 





The value of was obtained by plotting. Near the bus 
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bars readings of ee were plotted versus perpendicular distance 
from the bus bars. The slope of this straight line was see 
Bee Fig. 32. 

The determination of the velocity at a point on the blade 
required the measurement of aE at that point. A sample of 
recorded readings of ie is contained in Table VIII. These 


readings were converted to as 


ab 2 ra =o wien al 





by means of the relation 


dis the separation of the double probe. Values of al are 


also shown in Table VIII. Values of Ne were then obtained 


Vi 


by means of the formula 
= ab. coe |, taht randin i aa| can 

The NCR 102A electronic, digital computer of the U. S. Na- 
val Postgraduate School was used to compute a nel ; of, : of; ; 
and Ve . Machine assistance proved indispensible for the ex- 
Beive calculations required by the many cascade combinations 
tested during this investigation. 

One program on the computer was used to determine Cae Lins 
put data consisted of twenty potential readings, ten from the 
upper surface on the airfoil near the trailing edge and ten from 
the lower surface near the trailing edge, together with their 
respective distances from the trailing edge and the value of p. 
see Appendix E. The computer solved for on 


A second computer program was used to determine &, , oQ > 


and X3 - Input data consisted of er J oe : c. , Oe , © », ane 


(-&) Mlimeadition, inputs of © and ae were required for 
Qo 


35 





; : 
each of the points at which velocity ratios Ve. were desired. 
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7.  DIsSCuUssren 
7.1 GENERA 


The discussion of this method for predicting the flow of an 
ideal fluid across any arbitrary cascade first examines the re- 
sults of the tests made on various cascades. The analysis of 
the results then provides the basis for a critique of the method 
and for suggested improvements in equipment and in technique. 
Possible applications of this method to other aspects of engi- 


neering are also suggested. 


(PR Wasluseies 


The geometric parameters which determine a cascade are the 
blade profile, the stagger angle, and the solidity. For a given 
cascade the quantities K and Pr, are the fundamental parameters 
which were determined by this investigation for indicating the 
performance of a cascade over a wide range of mean velocity di- 
rections. Turning angles were also calculated for each cascade. 
From this information other parameters: were calculated to 
permit cascade analysis from different viewpoints and to permit 
comparison with published reports. 

It first seemed appropriate to compare the results of the 
method herein discussed with a published analytical solution. 
Ref. 7 contains an analytical solution of potential flow across 
a given cascade of { = 1 and 4 = 45, This report was used for 
comparative purposes because it contained coordinates for the 


profile used and because results were reported over a range of 
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directions of the mean velocity. The blade profile used is sim- 
jlar in camber and thickness to those presently found in compres- 
sors. The blade, however, has a trailing edge of zero thickness. 
The pointed trailing edge is difficult to duplicate exactly in 
the paper cutouts and would not be used in actual compressors. 
Reference to Fig. 21 indicates that the method of this paper ap- 
parently predicts coefficients of lift which are about 14% great- 
er than those contained in Ref. 7. The difference represents a- 
bout 5° in turning angle. This difference wae greater than had 
been expected, but it is subsequently shown that this disparity 
can be reduced. Subsequent tests indicated that an increased 
error in turning angle resulted from tests on cascades with large 
Stagger angles. The larger stagger angle cascades are more sus- 
ceptible to errors caused by the variation of specific resistance 
of the paper between width and length directions. Possible solu- 
tions to this problem are discussed later. 

Fig. 21 does indicate that the slopes of the curves of 
Kp vs.X are almost the same. Calculations reveal that the 
cascade of Ref. 7 has a coefficient K of .74 while the cascade 
coefficient determined by this method is .72. Thus if determi- 
nation of the angle of zero lift were made separately, the turn- 
ing angle predicted by this method would be within a fraction of 
a degree of the turning angles reported in Ref. 7. 

Fig. 21 indicates that increase of the solidity of this 
cascade to 1.5 changes the cascade coefficient to .49. The 1.5 
solidity cascade is not analyzed in Ref. 7. | 

Ref. 9 is a report of systematic two dimensional cascade 


tests of NACA 65 series compressor blades. Cascades considered 
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to be most typical of present compressor design were selected 
from the report for investigation by this method. 

The NACA test data was taken by varying cascade stagger an- 
gle. This was required because of the difficulty of varying the 
direction of the approaching velocity in the two dimensional 
wind tunnel. The method of this report permits variation of the 
direction of approaching velocity for constant cascade geometry. 
Comparison, then, was possible only at certain points, indicated 
Mme gs., 22, 23, 295 and 26. 

It must be borne in mind that Ref. 9 is a report of tests 
performed using an actual fluid. The magnitude of the effect of 
viscosity is uncertain, but one could not expect results based 
on viscous fluid flow to be identical to those for ideal fluid 
flow. Comparison serves to indicate, nevertheless, that pre- 
dictions based upon ideal fluid flow are sufficiently close to 
measured values for air to be of value in actual cascade design. 

Figs. 22 and 25 give the results for the cascade of the 
NACA 65-(15)10 airfoils with @ = 1 and 1.5 respectively. The 
figures show that the relationship between turning angle and cas- 
cade entrance angle is almost a linear one, indicating that the 
exit angle remains nearly constant. For a given cascade en- 
meance angle, turning angle is of course increased by decreasing 
the stagger angle. The effect of increasing solidity on turning 
angle (and therefore on the quantity TX) } ag small. Hence the 
cross plots indicate an appreciable reduction in (Cp with 
increase in solidity. These are all general results which one 
would expect. 


The cascade designer wonders what represents optimum condi- 
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tions and the extent of variation permissible for satisfactory 
performance. The cross plotted values of Ly give an tndiveas 
tion of the coefficient of lift associated with a certain turn- 
ing angle. This information together with the plot of ideal an- 
gle of attack helps to denote an area of satisfactory operation 
based upon upstream conditions. However, for the decelerating 
flow experienced in a compressor, information showing the extent 
of the range of operation for an optimum drag-lift ratio is de- 
sired. For this reason Figs. 24 and 27 have been presented. 
Ref. ll says tests show that the drag-lift ratios of cascades 
varying from turbine to compressor types show a minimum at an 
approximately constant value of cy, of about .8. For blades 
of the types generally used in cascades the drag coefficient is 
low and almost constant over a range which extends perhaps five 
degrees above and fifteen degrees below the minimum drag-lift 
ratio. These general criteria for determining range of opera- 
Groneare indicated on Figs. 22 and 25. 

The effect of changing cascade stagger is apparent from 
Figs. 22 and 25. One would expect the spacing between lines of 
constant 0 to be more uniform than this test indicated. This 
Variation is considered to be a manifestation of the experimen- 
tal error of the method. For the whole series of tests the error 
is believed to be + 1°. A test program which included a wider 
range of stagger angles or which provided for smaller test in- 
Crements in this same range would serve to assist in reducing 
this random error by permitting the application of statistical 
data smoothing techniques. 


There is again an indication, particularly for the cascade 
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Best at § = 1.5 and 0 = 35° that the high stagger ancuemeca= 
cade is subject to an additional error. The turning angles pre- 
dicted for this case seem to be greater than one would expect, 
considering the other three cascades. As in the 45° stagger an- 
gle cascade first discussed, some of this difference is probably 
attributable to the paper. 

This particular cascade was also analyzed to see what value 
of On would give results in agreement with the NACA report. 
The point on the blade cutout which had the value of je so de- 
termined was located on the upper side of the blade profile very 
near the trailing edge. This suggests that the real fluid in 
the NACA tests, under the influence of viscosity and adverse 
pressure gradient, separated from the airfoil at a point closely 
approximating this same location on the blade cutout. 

For the J = 1 cascade the predicted turning angles com- 
pare remarkably well with the NACA test points for a cascade en- 
trance angle XK, = 45. The comparison is not as close for the 
mee NACA points plotted at a cascade approach angle od, = 30%. 
hor © = 20° a line drawn between the two NACA test points at 

ol; = 30° and X, = 45° would not have the same slope as the 
corresponding line produced by this method. The same thing is 
true to a lesser extent for 2 =—257 seeimesemoatterences mum 
Slope are considered to be a consequence of viscosity. 

For the (Q = 1.5 cascade there is a somewhat larger differ- 
ence between predicted turning angles and reported NACA test re- 
Bets, but the comparison is still close. Twelve inch chord 
blade templates were used for three of the four cascades of 


( = 1,5, Only the 6 = 25° cascade was cut out using a $1ix 


£ 
A 
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inch template. The somewhat increased error due to bus bar 
location for the larger template is not discernable; so use of 
the larger blade is probably warranted for high solidity cas- 
cades when the reduction in cut out time and in labor are con- 
Sidered. A Series of analyses for various values of p and 
such as is contained in Appendix F for Q = 1 and 6 = 45° 
would permit approximate determine of the largest Sized 
blade permissable for a given acceptable percentage error. 

Fig. 28 presents a concise summary of the effects of solid- 
ity and stagger angle upon cascade performance. The variation 
of cascade coefficient K with changes in stagger angle are 
Shown for each solidity. For purposes of approximate comparison 
values of K for a flat plate cascade having the same nominal 
Meeue Of SOlidity and stagger as the actual cascade are also 
plotted and are seen to be GQuive “Similar. 

Another concise representation of cascade characteristics 
is contained in Fig. 29 which shows the effect oe changing stag- 
@emeangle upon the relative angle of zero lift o&,,. This pre- 
sentation amplifies and clearly shows the variation in the an- 
le of zero lift with changes in stagger angle. It suggests that 
a Similar plot for a much larger number of tests would permit the 
determination of a smooth curve of om vs.6. If the magnitude 
of random error were thus partially removed from the determina- 
tion of Za. » a more uniform spacing of the constant 6 lines 
in Figs. 22 and 25 would result and the error indicated in Fig. 
21 would be considerably reduced. Unfortunately time did not 
permit conducting the additional number of tests which would be 
required to determine this curve. 
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ye 
inescurve of Fig. 31, (| vS. percent chord compares favor- 
ably with the NACA curve for a similar cascade entrance angle. 
However, Fig. 3l is a smooth curve determined from calculations 


based upon raw data by means of a method of least squares de- 


- veloped by Professor Gawain. Before extensive tests of veloc- 


ity could be undertaken, the double probe and perhaps the tech- 
nique of ‘obtaining A= readings would require improvement. Fig. 
31 does indicate, however that the velocity distribution can be 


determined by this method. 
7.3 ACCURACY 


The accuracy of the results obtained by this method of 
cascade analysis is a function of the approximations made and 
Beene accuracy of the measured input data. The results are 
more sensitive to errors in some input quantities than to errors 
in others, and certain measured quantities can be determined 
more accurately than others. 

Table X gives an indication of the sensitivity of results 
to variation of input data. It lists the change in calculated 
turning angle caused by individually increasing each of the in- 
put quantities e, ce. Cas and 2G. by one percent of its 
measured value. This information was compiled for a cascade of 
the airfoil contained in Ref. 7 with 6 eo eand J = 1, bie 
it is considered to be representative of the other blades and 
cascades tested. The table indicates that the results are par- 
ticularly sensitive to changes of as and om. It indicates 
the desirability of obtaining the input measured values to an 


accuracy of one part in ten thousand. 
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Many operations were required for each cascade layout, and 
every step had some effect on the value of the final measured 
quantities. The largest of the construction errors was due to 
the conducting medium. The value of the specific resistance of 
this material was given by the manufacturer as 400C ohms. This 
value was speci.ied accurate from +3 to +8 percent with the 
lower value of resistance along the length of the roll and the 
higher value across the width of the roll. This variation is 
explained by the manner of manufacturing. These errors were 
necessarily accepted as the particular conducting medium used 
vas basic to this investigation. 

Wood was used for structural members of the working sur- 
face. The position of the bus bars was held to the tolerance 
mentioned in the equipment section, but the trueness of the 
joint depended upon the ability of the wood to maintain its di- 
mensions under continued use. It was found that after three or 
four tests, the table required cleaning because the graphite 
depo€éit on the table caused a change in the conductivity of Sub- 
sequent papers when not removed. 

In the layout of the cascade on the paper several sources 
of probable error were encountered, such as dimensioning, con- 
Structing layout and reference lines and setting stagger angles. 
These were subject to the same limits of error encountered in 
engineering drawings A reading glass and a Peive lene eyepiece 
were used to provide magnification for greater accuracy in lay- 
out. Dimensions could be held to a tolerance of .O1 inch. 
Lines were maintained to widths of from four to six thousandths 
of an inch and angles were maintained correct to 10 minutes of 
rc . 
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The templates introduced a systematic error due to the 
tolerance of manufacture and a random error in the process of 
positioning. The manufacturing tolerance was + 0.005 inches. 
The profile cutouts were maintained-to their respective coordi- 
nates with a tolerance of about .O10 inches. 

Reading and determination errors of the instruments were 
encountered only with the null meter and with the probe. Since 
a zero reading or null was required on the meter prior to read- 
ing the voltage divider, detecting a null was reduced to detect- 
ing motion of the meter needle. The sensitivity of the instru- 
ments permitted reading potential on the voltage divider to one 
part in ten thousand. Readout was directly in percent poten- 
tial; so constancy of imposed voltage was not required. Errors 
in probe location depended upon all previous dimensioning and 
construction of reference lines. With optical aids it was pos- 
Sible to place the probe within the actual width of the pencil 
lines. Distance from the cutout edges of the profile could not 
be determined to comparable accuracy. However, the constant 
potential lines are normal to the cutout near the profile edge; 
Memtnis distance is not as critical. 

The velocity distribution around the profile was determined 
to demonstrate in preliminary fashion one of the capabilities of 
the method. This involved measuring small changes in potential 
and distance as a means of determining potential gradient. 

Tests were not conducted to determine the accuracy of these 
measurements nor the validity of the final results. Ref. 12 
States that for a liquid electrolyte it is possible to obtain 


an accuracy of two to five parts in a thousand with respect to 
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the theoretical value when using a double probe of five milli- 
meters separation. The method of determining the error of such 
an instrument was to rotate it around a point in the field in 
which case a sinusoidal variation should result. 

It might seem, when first considered, that a comparison of 
the potential readings at corresponding points on adjacent pro- 
files would give an indication of the overall accuracy actually 
obtained. However, this is not a definite criterion. In fact 
it was not uncommon for the measured potentials at correspond- 
ing points on adjacent profiles to differ by 2% without produc- 
ing a noticeable change in results. The reason for this is that 
the quantities TC and ‘ which determine conditions in the pic- 
ture plane contain ratios of differences between measured quan- 
tities. The measured values of the various €'s could there- 
fore differ from profile to profile without necessarily affect- 
ing the values of CT and £ 

A better indication of the overall accuracy actually ob- 
tained in construction and measurement was the repeatability of 
Beewets. Fig. 33 gives an indication of the degree of rerpeat- 
ability of the method and procedure as finally evolved. The 
three layouts shown were distinct in every respect even to the 
use of Teledeltos paper from three separate rolls. 

The indication of repeatability does not, however, purport 
to show the accuracy of the method in its degree of approxima- 
tion to the true conditions for an infinite cascade. The bus 
bars were located a finite distance from the cascade and can on- 
ly approximate the conditions an infinite distance away. Al- 


though the error accepted with this assumption cannot be deter- 
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mined exactly with the present equipment, a close approximation 
of this error based on a cascade of flat plates is derived ana- 
lytically in Appendix F. From the figure and the data in this 
appendix it is evident that for the geometry selected, the er- 
ror introduced by the location of the bus bars is negligible. 
Only one cascade of flat plates was investigated, but the ana- 
lytic means is at hand for determining cases of different 
solidity and stagger angle. 

The effect of the foregoing limitations on the accuracy of 
K and ©. cannot presently be determined directly. Further 
testing is needed to provide statistical data for determination 
of the error as well as means for its reduction. 

The closely related parameter of cascade performance GO, 


when compared to published data, agreed within the limits of +1°. 


7.4 CRITIQUE OF THE METHOD AND SUGGESTED IMPROVEMENTS 


The given conducting paper was taken as fundamental to this 
investigation. The scope of this work could not include a 
search for nor an evaluation of various conducting mediums. 
However, removal of the directional properties of the conducting 
medium and reduction of the variation in its specific resistance 
Should greatly improve the accuracy of the method. Perhaps the 
present type of conducting paper could be improved. Other con- 
ducting mediums may be available. Considerable work has been 
done in France and at Rensselaer Folytechnic Institute in this 
country with electrolytic tanks. See Refs. 12 and 13. A sol- 
id conducting medium is believed to be preferable for the meth- 


Od herein described. The solid medium obviates such problems 
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as leveling and polarization. The test set-up is much simpler 
as only a single template is required oe all of the cascades 
Meyng that airfoil. 

If after investigation the teledeltos paper still proved 
to be the best solid medium available, allowance for the di- 
rectional properties could be made by use of a scale factor. 
For the work reported here a scale factor correction of profile 
coordinates was not considered practical because of the range 
of variation of specific resistance in a given direction. Sum- 
marizing, even if directionality cannot be removed, variation 
of specific resistance in a given direction should be decreased. 

A minor suggestion for increased ease in cascade layout is 
to change the black color of the backing side of the sheet to 
more clearly show-up the penciled layout lines. 

The layout table undoubtedly would be further modified if 
this method were to be used commercially. It is considered de- 
sirable to bring the meridional ends of the paper together to 
establish an electrically continuous sheet although the origi- 
nal flat table produced satisfactory results when the approxi- 
mate streamlines were used as boundaries. The cascade test at 

(7 = 1.0 and © = 30° was done on the original table. If a 
paper vere ee which was sufficiently long in the direction 

of the cascade axis the effect of straight-cut ends would prob- 
ably be small at the center blade cutouts. This would be par- 
ticularly true for the cascades of small stagger angle. 

Miva table be constructed of design Similar to the one 
herein described, particular care should be taken to insure 


that the table hinges be dimensionally true. Their trueness 
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also will require periodic checking because of the dependence of 
the electrical connection upon them. 

Perhaps a clip-board type of springed clamp could be design- 
ed for establishing electrical contact between bus bars and pa- 
per. Such an arrangement could save time over the bolted clamp- 
ing system used for these tests. 

Another suggestion for improvement is to design the layout 
table such that the bus bars as well as the paper would be elec- 
Memedily continuous. Although the conductivity of both the cop- 
per bus bars and the silver painted areas on the paper is large 
relative to the paper, their conductance would improve the ap- 
proximation of uniform conditions at large distances from the 
cascade. 

The method is dependent upon the technique of the operator. 
With an increasing number of cascades tested, the operator can 
be expected to improve in accuracy and increase in speed. The 
angle of holding the probe, the depth of probe penetration, the 
eyesight of the operator, and his ability to position the probe 
accurately are all personal factors which will affect the re- 
sults to a certain small degree. 

The statistical curve fitting technique for determining 
the value of On probably could be improved. There is a re- 
quirement for reduction of the random error. The method used 
in this report seemed attractive because of the similarity of 
its equations to those for perfect fluid flow. It proved to be 
practical for solution. This same form perhaps could be modi- 
fied to allow for the curvature of the blade profile near the 


trailing edge. 
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This method of cascade analysis has the great advantage of 
being much less expensive than actual cascade tests. The 
teledeltos paper is the principal consumable item. The templates 
can be manufactured without difficulty or excessive expense. 

The time saving of this method over other methods of cas-~ 
cade analysis is considered to be one of its principal advan- 
tages. With improvement in technique it was found that a typical 
cascade layout could be completed and measurements taken for de- 
termination of turning angle in slightly more than one working 
day. This of course depended to some extent upon cascade geom- 
etry. Calculations of turning angles from the raw data were ac- 
complished in three to four hours and would be much faster on nev- 
er digital computers. The use of the computer is not required 
for this method. It was used here principally to save time and 
thus permit the testing of a greater number of cascades. The 
measurements required for the determination of velocities at 60 
Stations on the profile required an additional two hours. Cal- 
culating time to find these velocity ratios required approximate- 
ly 20 minutes for each angle of mean velocity. 

Analyzing the method from another viewpoint, it experimen- 
tally accomplishes conformal transformation of an arbitrary cas- 
cade onto a circle. When one considers the well known transfor- 
mation vhich maps a circle into a cascade of flat plates, he 
might in one sense consider that this method determines an e- 
quivalent flat plate cascade for any arbitrary cascade. The 
performance of a cascade of flat plates is defined theoretically 
by two parameters K and Ol . These quantities are determined 
in this method simply by making the measurements required to 


mind #3 and fo 
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This method has greater capability, however, than simply 
the determination of these two parameters. A point to point 
correspondence between blade profile and circle is obtained 
which permits determination of the local velocity ratio at any 
point on the blade profile. 

In review, a method of cascade analysis has been developed. 
It was shown to be capable of predicting, with reasonable accur- 
Segeche turning angles of perfect fluid flow across an arbi- 
trary cascade. Because of the nature of the present conducting 
medium the accuracy seems somewhat poorer at large stagger an- 
gles but no limits of solidity or camber are known to exist. 
The time required to accomplish the procedure is short enough 
to make the method suitable for engineering application. The 
method is not expensive. 

Extension of the method to the analysis of radial and con- 
ical (mixed flow) Paecadee is analytical once the basic trans- 
formation of a rectilinear cascade has been achieved. 

The method could be analyticaily supplemented to permit 
consideration of viscous conditions. Some estimate of com- 
pressibility effects could be made using the Prandtl-Glauert 
method. | 

The method of experimental conformal transformation used 
herein should be a valuable tool when applied to the many other 
engineering problems which are characterized by LaPlace's equa- 
tion. Examples are problems in the theory of elasticity, in 


heat transfer, and in electric field theory. 
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8. CONCLUSIONS AND RECOMMENDATIONS 


The principal conclusions of this investigation are: 


The method herein described determines the ideal fluid 
flow, including turning angles and pressure distribution 
along the blades, for any two dimensional cascade. A 
practical technique, partly experimental and partly ana- 
lytical, of conformally transforming an arbitrary cascade 
to a corresponding circle is provided as the basis of the 
method. 

At its present state of development the method appears 
to determine turning angles to an overall accuracy of +1°, 
but further refinement can undoubtedly reduce this error. 

The cascade analysis is accomplished in approximately 
one working day and at nominal expense. 

No limitations in the method due to cascade geometry 


are known to exist. 


For further refinement of the method it is recommend that: 


A more nearly homogeneous and isotropic conducting med- 
ium be used. 


Cascade layout and measurement procedures be mechanized. 
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Potential supply 


Null meter 


Voltage divider 


Potentiometer 


Voltmeter 


TABLE 
INSTRUMENTS 


Heathkit Battery Eliminator, Model BE-4. 


General Electric DC Micro Ammeter, Type 
DO-71, Model A29AA1, coil resistance of 
1700 ohms, detecting sensitivity power 
on the order of 8 micro-micro watts. 


Electro-Measurements Inc. Dekavider, Mod- 
el RV622, linearity of +0.0025 percent, 
input resistance of 10,CCO ohms, Serial 
ba, 04ce 


Leeds and Northrup Portable Precision 
Potentiometer, Number 8662, Serial 
13LS7662 . 


Weston DC Voltmeter, Model 45, Serial 
49400. 
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AIRFOIL ORDINATES FROM NACA ARR. NO. L4K22b 


TABLE- ia 


(Stations and ordinates in percent of airfoil chord) 
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TABLE III 
ORDINATES FOR NACA 65-010 BASIC THICKNESS FORMS 


(Stations and ordinates in percent of chord) 


Stato, xX Ordinates, x+y 


65(216)-010 airfoil combined 
with y = 0.0015x 
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TABLE 1V 


COORDINATES OF 6" and 12' NACA ARR NO. L4K22b AIRFOILS 


6" Chord (Stations and onrdimapes in inches } 
Upper surface Lower surface 
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12" Cherd (Stations and ordwmaccemiunmeumenes } 


TABLE 1V (Cont. ) 


Upper surface 
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9.6348 
Poe 22 76 
16.8204 
4 120 
m2. 0000 


x Ordinate, 
0 


a oe el en ed en oe eed en 


Ae Bey: 
. 1488 
5 EIS 2 
Hoy aye. 
9004 
Tone ee 
ee 
roo 
1976 
348 
»-4640 
.0460 
-6068 
mere ee, 
.6380 
6044 
Raye Cie, 
»4304 
Pe 
. 1436 
9648 
» 1680 
.0472 
. 3084 


Q 


Y 
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ee 
NFR OOOMOANYNAOOD SWWN RF 


Power surface 


Diver one x 


Q 


0924 
212 i2 
» 1944 
3564 
6684 
.9756 
2192 
so 192 
.4744 
Helo mir 
»6540 
2432 
seele 
4216 
. 0048 
Se lle 
LG 
A levers 
. 3640 
276.16 
Oo. 
1724 
» 1796 
» 3550 
0000 


Ordinate, 
O 


rae 2 
pOOZG 
BODO 2 
0444 
301 D0 
FO ZES 
OTs, 
» 1296 
 l96e 
mee 
28 
soe 
» 3948 
»4296 
»-4692 
PoOGe 
S5oee 
palo} ic 
so i2a 
5604 
oe oe 
4728 
P3720 
RP LPA ONS, 
Q 
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TABLE V 
COORDINATES OF 6" and 12" NACA 65-(15)10 AIRFOILS 
6" Chord (Stations and Oneum-scomnmeinmehes } 


Upper surface Lower surface 
Seat lO met x Ordinate, y Steimon, xX Ordinate, y 
0 OQ 


0059 0606 s0542 = ,O FSC 
oS 1 207 10 .O719 -.0146 
.0438 . 1079 Coz eo OlG 
mel 22 ley OL . 1878 Be O2 4 
Saaye Pose » 3442 +,0164 
»-4033 » 3465 4967 Used 
5 Se pse .4132 .6474 Oa eal 
.6544 a aa, 9456 »0528 
ie LOS 7 .6079 Laat » L055 
1.4643 .6746 BP aye rsle. mleo? 
wo 7/09 pio loco! 1494 
2.0780 Oa Pole eo LOO 4 
ooo S peaere)® 2.4147 . 1789 
2.6934 » 1944 aot 2. ho tal 
3. J000 . 1886 370000 »2042 
E2066 . 7/009 S27 a4 P2166 
0120 7 307 Selec toe oe 
e27160 Be OS 3.8640 »2474 
4.2184 sOnlee 4.1616 2G 
4.9191 0456 4.4509 Sale) 
4.8181 »4675 4.7819 » 2489 
ee 152 Bowie): 2.0546 »2294 
9.4114 2/0) Ore ielele » 1892 
e. 70/0 O22 See, sree 
6.0000 . OO90 6.Q0UG0 = 0090 


LEADING EDGE RADIUS 


= ,03996 inches 
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be 
OOOO MANANABOOHBDHWHDNMEH 


b+ 
NO KF 


12" Chord (Stations and ordimacesma aeenches } 


Upper surface 


Station, 


0 


POLLS 
I SO2 
BOG /6 
»2245 
mo 116 
. 8066 
BLOODS 
BOSS 
3174 
22207 
pod l/ 
wD D9 
er 7O/ 
. 3868 
0000 
pOlS 1 
»2241 
Boo!) 
-4367 
BUSS2 
oycleyP2 
» 2304 
pe2 28 
414] 
.0000 


X 


Ordinate, 
O 


ed ed ee ee oe ee ee ee eo 


TABLE V (Cont. ) 


ees 
Miele) 
a eaks S 
340] 
omy 
»6930 
.6264 
.0459 
Pe ioS 
3492 
~4507 
FOZCS 
7271 00 
POG SS 
ed al 
eee [ 
»4614 
sol 92 
pcos 
nO 
A erS Ie) IL 
iD Zan 
some 
» 3249 
20180 


ye 


BEAOLNG EDGE RADIUS = 


be be | 
DOOMDMANANAOOHBDHAPWWHrH 


ee od 
NO Fr 


Lower surface 
Ordinate, 
O 


Slawron., xX 
C 


OeeZ 
~ 1438 
.2124 
wo 1oo 
.6854 
9934 
.2947 
ended 
4826 
OS 
Gays) oie! 
~2440 
oes 
20344 
0000 
PP aielexs! 
~ 1759 
Bory 
ROO OS 
5 Shel les: 
yeejele 
. 1696 
eae 
See) 
»O0O00 


.07992 inches 
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Soret FE 


Osis 
O2Z72 
O253 
0053 
0328 


NOWAG 
» 1042 
OOo 
ae ha 
i2ole 
eee 
3334 
MSS T7s 
Bley 2 
.4063 
-4373 
.4664 
-4949 
Oe 
ap ee 
4977 
.4567 
~3784 
2443 
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TABLES 
TEST PROGRAM 


fhuching: Amoles ; 
l. First series: 
Profile - Mutterperl (NACA ARR No. L4K22b) 


R 
Tests (a) © = 450.) 10 Seaeee 
[b) 6 = 459. 7G) ieee 


ee Second series: 


Profile - NACA eee 
Tests - Group one (six inch chord) 


(a) 2 = 200 Gale © 

(ob)? = 250. 0 

bs © = 300-5 Gf = 1.6 

d) 6 = 359. sie= io 

Group two : 

ha) © =" 002 lo ie inch chord) 
to) 6 = Oot mie om ci) tnehmenord) 

‘¢ 6. = 2 3C2as Ue Hl Sane he inch chord 
d) © = 220-= 0 =sles (twelve inch chord 


Velocity Distribution: 


Profile - NACA 65=(15)10 (twelve inch chord) 
Test - 6 35°: § = 1.5 
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TABLE VII 
SAMPLE DATA FOR TURNING ANGLE DETERMINATION 


Upper Lower 
Surface Surface 
Station € (Percent) Station € (percent) 
3b 71.88 46 30.69 
S Fee sts: a 30.6u 
d 71.89 b 30,85 
4 71.90 C 30.49 
4a 71.89 d 80.42 
b eos 47 ClOm ei) 
a 30.29 
b 3C24 
46 34.91 Cc 3073 
a 34,70 d 30,412 
b 62.47 48 30.08 
c Ba..26 a 30-03 
d 64.00 b 305Cu 
47 Sani 1 € 29.99 
a 6.52 d 30.00 
b Gen23 49 30.18 
E Selon 
d 32.76 (3&) 
48 32.50 Kilo * - 
a 82,20 
b 31.93 Station § C(Percent) 
C Si O9 
qd ei, 24 O 98.53 
49 S663 1 95.74 
Center Tail 30.52 2 92.97 
S 90.24 
4 G1 oe 
s 84.84 
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TABLE VIII 
SAMPLE DATA AND INPUT VALUES FOR DETERMINATION OF VELOCITY RAT10S 


Upper Surface Lower Surface 


ob 
Percent CS nis ay e aE Ag 
Chord (Percent) (Woelts) (Percent) (Percent) (Volts) (Percent }) 
meom L.E. 
O Be cc 67 a64 NAN 

2,05 fale, Sl LOR Zz 2235 S30) C4 .cG 7.963 
on 1 / (Lt ga) S206 0.824 63.67 29 288 6.732 
ee 2) Wess - = 62522 De ES) 6,312 

E33 (Oye Zao 0.592 fe RO) @ 22207) 5.266 
mo. 42 mn 53 vA ae 1.098 59.85 27 Eg A, 4.995 
2. 20 Wale. 23 6.06 1.395 58./4 19.71 4.5937 
i. 67 70.40 8.91 POS |. slopes ll 19.29 4.441] 
2O.33 69.29 10.69 2746 ll 54.05 \RsPar ee) Al SNe 
20,00 67.99 Deke S020 O2 om 18.70 4305 
29.17 Bor OL 13.93 Src O7 51206 16.30 Shoe 
oo. 33 64.88 Loi 6 Soe 49.26 LOGS Bee 
ey. 50 eo.06 Lees 4.197 Na Silale, bem ss 35026 
41.67 one, 2 1 18.43 4.243 45.95 16.225 3.741 
7.03 59.18 VS eer Aes 2 44.34 ey Sees 
boO.00 ayers ne. 20.31 Ae tee: 42.74 Nina wa 3.630 
a. 17 55 06 20.12 4.632 Ano 14 7 BSe4 
mo. 33 Bes IS 19.64 4.521 39.81 see he ls S3027 
62.50 70.6 / 19.53 4.496 38.44 13.08 Cruep ml 
66.67 ao oO) Teese L213 Cee 123 2cOZ 
fe .03 AG. OG 18.98 4.369 35.94 VAs 2.090 
fo .00 44.57 19.40 4.466 34.79 elo! 20am 
79.17 A OZ e349 Ane 7 Soe 2 Loe OG 22302 
B3.33 40.39 13234 A222 S2e 7 9.22 Oe 
B7?.20 So.23 19.3! 4.445 3169 6.21 1.890 
31.67 Se. 10 20)..09 Ay Feil B1.05 Sees’ 12565 
93.75 34.91 20.02 4.609 30.69 5.94 l2ser 
ao. OS 33.97 19.85 AoW Ae) 80257 oo He) L232 
97.92 S250 20470 AGO 20,05 Cee) O2eS4 

Az = .0405 aches 


OM 








bb 


TABLE 1X 
INPUT VALUES FOR DETERMINATION OF TURNING ANGLES 


Ne, 
im & a Ga ee Pia s 
Deepees —-Lercent rereene Bee Ot jevalelayeae LDches 
MUTTERPERL AIRFOIL 
0 45 61.86 39.42 40.29 2.870 6.00 
Be) 45 S261 38.02 38.54 2.825 4.00 
Ag 20 Sulpateke: 40.54 40.63 3.060 6.00 
0 25 el.14 . 40.32 40.50 32072 6.00 
O 30 62.14 39.01 39.30 S23 259 6.00 
0 35 62.34 40.82 41.26 Lee 6. OG 
5 20 . 71.38 S0n26 S066 2.943 8.00 
i 25 62,63 39.70 Bones 3.019 4.00 
5 30 fe. 26 30.03 30.44 Be Tec 8.00 
5 35 717.90 29.99 30.49 LPT eal 8.00 


el ll ee ed 
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TABLE X 
BENSITiIVITY OF RESULTS TO VARIAIR@eiS Or INPUT Valves 


Quantity Value MQuantity OQuantity at pts 35° 
Qua hie Ty 
6; 62.91 +6 001 + 012° 
oF 38 .58 +4 .001 ph Melo 
oe 39.68 +4 001 aie loac 
2 30.80 +3 001 + 029° 
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atom= 60° 
+ .026° 
+.641° 
~.647° 
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Figure 5 -Simple Two Dimensional Flows 
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Figure 6 -Correspondence Between The Real and Picture Planes 
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Figure 20,— Cutting knife modifications. 
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APPENDIX A 
THE ANALOGY 


In potential, two dimensional, ideal fluid flow there is a 
complex velocity potential w = téy such that the Cauchy- 


Riemann differential equations are satisfied 
Pp . 2dY . 7 § 
and 
P —= _ sh 
aN og 


Then the complex velocity may be written 


E 5 oP. : : 
He P-GE edpewie 


Considering the flow of a steady current in a flat sheet of uni- 


Al 


form resistance per unit square V . Y 


. va 
Cn =a % A3 
Pa yn 
where én is the current line intensity. Let 1 be a total cur- 
rent function analagous on the sheet to the stream function in 


the ideal fluid alton such that 


= Nel ; ue 
Re oy a : 
ice tOollows: 
Vi 2 yl _ ME 
lar smecen and wv Sy ; Ad 


which indicates that E and I also satisfy the Cauchy-Riemann dif- 


ferential equations. Equations AS may be used to obtain 


E +épl Cal) = 2), Me 
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or if we divide by the maximum potential 





F* = Re AT 
E+ iti eCrid = oa Ag 
Hence we may list the following analogies: 
Z Plane ¥ Plane 

Meee = 82 . £* fae Ba 

cm v ox dy v Cs) oe AQ 
2 - 1 aL. £*(_ 26) yp Sele Senora 
mer dy) ox PN By o7.COE 
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APPENDI AE 


DERIVATION OF EQUATIONS FOR CONFORMAL TRANSFORMATION 


This appendix develops the equations needed to achieve con- 
formal transformation of the cascade of blades in the z plane to 


d 


the circle in the 5 plane. It also defines the quantity di 
which is important in later work for relating velocities in the 
picture plane to those in the real plane. 

Conformal transformation is possible between the real, z 
plane and the picture, 5 plane if there exists a function in 
the picture plane of the real plane complex variable z such 
that for every assigned value of the variable the function has a 
definite value or system of values and also a definite derivative. 
The advantage of this definition is that it does not require the 
existence of an analytical expression of the function. 

lt was shown in the main text that any point on the real 


plane may be designated by its coordinates x, y or by the values 


of o and wt which characterize it. See also Appendix A. Fur- 

:. bee 
S \°3t/00 

invariant quantity which must be maintained to achieve conformal 


thur, the parameter was shown to be the unique and 


jmanstormation. 


The basic potential relations for a point of strength m are 
Bey gat =- 2 
p ee Ink ond YO Bl 


where mis real. The point is called a source if m is positive 
and a sink if m is negative. R is the radial distance from the 
point and Q is the angular coordinate with respect to the 


Bernt aS Origin. 


iGse 





In the picture plane of Fig. 4, let sources of strength m 
be located on the real axis at E=—b,-O%, and let sinks of 
equal strength be located on the real axis at&=b, a, 


Then the complex potential in the picture plane is 





m jy bt3 sy +p _ 
ur = = In Doe ins ~ oY, = (Oc. 0 Ue B2 


Now any point in the picture plane can be designated by its coord- 
inates § »Q or by the values of vy and WY which characterize it. 


Combining the factors of eq. B2, 





wr =D (np bt $4 , B3 
é7  b-3 ‘3-% 
On the circle 5) = 
wre = In ae’? fae? 4 ae’? ++ %V% 3; 
ra -qaeiS Facts + ae -t9 oe % + op) 
and with be = e® 
Ween co cosh¥+acos@ | 
in cosh¥ -Qcos@ BS 
We may let a = l, giving 
ur =F tn she + COs - o+iy B6 
at cosh¥ -cos@ ~ 
which shows that the circle is a line of constant W . 
If 1 and 2 are the points indicated in Fig. 4 
D-D = ncosme te | 87 
! ” Coshy — | 


By Appendix A, Pp corresponds to CS ; so we may impose the rela- 


tion 


p=k,-—k,© ’ B8 
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where kK, and ae are arbitrary constants. It folloveeenae 


P.-P, = K2 (eee) B9 


and if ies sk; (-3€ | ’ BLO 


i 6) 


P.-F _ uk €,-€,) : 
then _— ae =e bi = SE Bll 


This establishes the existance of a function relating the z and 
5 planes conformally. 


Eq. B/7 may be solved to give 


= 7 _| ; 
¥ = cosh™ € ae =. ee eines SLE) 
ff ett _ | Ble 
e — | 
Knowing @ , the locations of the sources and sinks in the pic- 
ture plane are determined. For any point P on the circle and a 
corresponding point P' on a blade profile, let there be a quan- 


tity £ such that 


i a ee B13 


PE, -f Cae, 
By substitution of eqs. B6 and B7 into eq. B13 we obtain 
in coshy + cos in coshat/ 
fod sh¥-cosG_ t*" cosht=l ae 
KR in cosh® + | 
COSNG = | 
and solving for cos@: 
f=| 
(cosh!) _ at (2£ - 1) | 
Bes = coche |e = cosh = B15 


(soso | et tAS~!) | | 
ie 


itesimal vectors in the picture plane to corresponding vectors 


The derivative is the operator which relates infin- 








in the real plane, and vice versa. In particular, 


Or 





Vp = gu = ache = =e = B16 


but recalling our parameter ‘ this may be written, 
V2 = Vala dz) B17 


Nov. V& is known in the picture plane and is given by eq. C5 of 
Appendix oF 


4 


The numerator of eq. Bl? can be vritten as 








df} x - B18 
d2 aa: > ; 
The denominator of eq. B17 can be written as 
ef : sin cosh® - 
~ | (cosh”¥- cos* osh T+] 
| COS te eae 1 
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APPENDIX C 


DERIVATION OF FLOW EQUATIONS 


This appendix derives the equations required to develop an 
arbitrary flow about the circle |5 |= O. 
The sources and sinks used to accomplish transformation can 
be considered to create an axial flow velocity Ve such that 
= \Q8 . On the circle 1 | = Q this flow produces a ve- 


locity. From eq. B3 


z 2£- Ves - 2\es. cosh? sing 


TF cosh2et—-cos2zg °* Cl 


The basic velocity and stream potential relations for a point 


vortex of strength AN are 


p-L-¢ and sg lek : C2 


Cal 
where AN 1S positive if counterclockwise. 
mr vortices of strength J\. and with direction indicated in 
Fig. 7 be placed in the picture plane in addition to the sources 
and sinks the resulting complex potential is 


. Vas b+3 4+, nye b+ y 
UW Vas (i BE on) igp(n Bey D-¥ + ind) Gs 


and the corresponding velocity on the circle is 


__ 2\as _cosh%sing 2/\. _sinh{os© 
Ve = T cosh2r¢%-cos2@ ‘ 7T cosheae—-cose6 ~- o 
If Va = and Vi = ; 
Me = tan Xm = es CS 
Vo 


Since the value of m is determined in the process of transforma- 


Oo 





tion, from m = $k,(- xe , the arbitrary choice of a mean veloc- 
Selo 

ity direction olin determines the strength JN...) Then if ad- 

ditional vortices of strength Jf and with direction indicated 


in Fig. 7 be placed in the picture plane, the resulting velocity 


Ji 


2 2Masfashtsing = —tandm sinh¥cos@| [7 Z%. 
a. Ve = coshet- cos 5 -E las COs 


At the point on the circle where O= On Ve =O > solving tog 


on the circle Sel. 1S. 








oi i COoSQw _ sinOr 
2m ’ Bs = (ON dm one ~ sinh F A 


The circulation given by this equation insures satisfaction of the 
Meeea condition. 


The velocity triangle is now complete. 





tan G = tan&m rae C8 
tan og = tand&m in] C9 
t a A * ; a 
an(ei-da) = eT a 
= tan@. 


Of aerodynamic interest is the ideal angle of attack. At this 
angle the front stagnation point is at the leading edge of the 


airfoil. IfV=0 forOG=O, equation C6 may be solved to give, 


. _ sin@r- sin€ ; 
ton ci = FanhG(coSGn ~ COS, a 


Fortine determination of lift coefficients, 


L= ty 4 Vr = PVnT . cl2 


From eqs. C8 and C9, a as(tan ol, — tana, ) Gis 
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Substituting and solving, 





Cp = 24 mlas (tana, — tano) C14 
_ XV | | ¢ 
put Vm = ree and m=VaS 3 $0 Ge Ala. ram $e oe (tanh ~ tano’) 
Hence, 
VO Z <08°K% (ran — tanoG ) cls 
Ly Gg cos am ee) 





APPENDIX D 


DEFINITION AND EXPLANATION OF THE TERM "SPECIFIC RESISTANCE" 


Consider an element of the conducting medium of length Al 
and width OW. Assume constant thickness h and resistivity f. 
Let AV be the potential drop along the length 41 of the ele- 
ment, and Ai be the cCOrrespondifigpeunrent. 


1+ 41-4] 


vt 






[a cele ae 
aa 


Ohm's Law may be written in the form: 

AVAL = He AI/ OA = (PAl(O%Aw) where AA= haw) . 

The constant of proportionality P/n is defined as the 
Specific PeeeGe atic é VY, of the conducting medium. Solving the 


above equation for this term: 


Y= (Mr) = AVAt (AMA1) wee tte = ons, 


amp. in 
The net resistance of the element may be written: 
A= AN /s? = v(4Yw) ; 
Thus if the element is square (Ql =A), its net resist- 
ance R, will be numerically equal to the specific resistance V, 


of the medium, regardless of the size of the square. 


eee 





APPENDIX E 
DETERMINATION OF Gr 


Lee we, =Az™ =6 ius e. 


be the complex potential for conditions symmetrical on the top 
| 


and bottom of a wedge, Fig. El, where Ry = i-z . EZ 
The sides of the wedge are two separate lines which join at the 
wedge vertex to form a single line bisecting the obtuse angle 


formed by the vertex. Along each line wWe=0, SO; 
x Re 
o Sala ee 


ie 
Similarly let Uy, = Bz” E4 
be the complex potential for conditions where a single line of 


constant SU coincides with the top and bottom boundaries of the 





=e 

wedge. Here pee = . Fy = F | BS 
aT A Fo 

_and as above, ea = +6r E6 


on the wedge profile. The + sign is for the upper side of the 
wedge and the - sign is for the lower side of the wedge. 


Superposition of these two complex potentials results in 


een Cc, Ar oe a 

Let the wedge approximate the trailing edge of the airfoil 
being investigated. Let r be the straight line distance of a 

point on the profile from the trailing edge. With the trailing 


edge at the point r = 0, eq. E/7 will result in € being zero at 


r = 0; so with the addition of a constant, the resultant equation 


is, oe = Cx a Alte | eee ES 


heise: 





Desired is the determination of the constants Cae A, B, 
and p, such that the curve of the above expression will as close- 
ly as possible approximate the measured values of © in the 
neighborhood of the trailing edge of the airfoil and permit de- 
termination of On The constant p is determined by a plotting 
procedure described later in the appendix. The determination of 
i is accomplished by the method of least squares which is 
mathematically expressed by making the sum of the squares of 


the deviations of the experimental points from the analytical 


/\ aa ie ~ é,) = minimum. E9 


curve a minimum: 


This requires that, 


DA. f -dsA-dA . | 
ye =SA YB _ 


The resulting set of equations is 
ay Ce tanh + ays B M, 
ae 7 + a»A bhp M, EDy 
aiCp +O2A +028 =Ms; 


) 


where 
: N N 
Q _ Q oe te le M = 4, Gs 
n=l, 22 ee i= ; 
éal é=! 
N N 
ae 2P - a _ | 2m 
Qi = 02) = 2M ’ aus? ga = 4 Den Me=W Lt a, E12 
=] = é=! 
N N 
7 mee! oF _ | 2 | AU 
Ai3= Os) = 24h", a= DN i M s=y dn nC, 
= e} ec} 


and N is the number of experimental points. 


ire Matrix solution for er is obtained, yielding 


c= M (0,,035— ais) +M —— ee +M 3(a i223 ~ Az2Qs E13 
= , 


ZO QigQe25 TF A2ds5 — 5 Qi2 C33 — OV C2 
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For this investigation a sketch of the trailing edge of the 
airfoil was drawn to a scale ten times that of the 6 inch chord 
templates used. Straight line distances were measured on this 
Sketch from the trailing edge of the airfoil to the locations a- 
long the top and bottom surface where potential readings were 
taken. A curve was plotted of en and G, versus these distances 
r: See Fig. E2. From this curve the values of Cy and C) were ob- 
tained at equal values of r. On log-log coordinates (Cy — Cy ) 
was plotted as a function of r. The slope of this straight line 


is the coefficient p. See Fig. E3. ‘ 


iS 
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APPENDIX. EF 


APPROXIMATION OF THE ERROR@EUESOsiheE BUS BAR Loetiia 


In adapting equipment to provide the electrical analogy to 
flow through cascades, it was necessary to locate the bus bars 
a finite distance from the cascade. Thus the location of the 
bus bars introduced an error in the measured data. The deter- 
mination of this error cannot readily be made when the airfoils 
of the cascade have camber and thickness. However, a good ap- 
proximation of the error may be found by considering a cascade 
of flat plates. 

Tests of this investigation were conducted under conditions 
equivalent to straight through flow with no circulation. The 
Same conditions may be achieved over a flat plate cascade by 
imposing a potential flow composed of sources and sinks onto the 


conformal map of the cascade. 


THES eASeApESOr FLA] PLATES 











Considering wy = 5) as a complex potential function, 
then 

c038 | 4 8tV%b sing as tn $2 Vb 
We ss 1 ager + 2S ing yao, al 


gives the mapping function between a circle with the radius a 
in the 4 plane and a staggered lattice composed of straight 
airfoils in the y plane (Fig. Fl). 

The first term of the above equation can be considered as 
originating from sources and sinks of equal strength located as 


Shown in the figure, while the second term originates from 
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vortices of equal magnitude superimposed on the sources and 
Sinks. The sense of the vortices is shown by the direction of 
the arrows in Fig. Fl. 

The spacing of the airfoils along the cascade is a con- 
Stant s and 6 if the stagger angle of the blades. 

The parameter b is connected with the chord Spacing ratio. 
In order to establish the relation between them, it is necessary 
mamcalculate the chord. By setting? =0 fon G=,iGeng.F !) and 
5=4r(T. Fig. F |) which points correspond to the leading and 


trailing edges of the blade, it can be shown that 


ke 
en. & (Yb + 45 cos” BP +))* + &y cos 
a cos® Ln e 


‘ (- / 


, -"| % sin 

+ sin® tan igs z/rcosB TIF 

For any chosen stagger angle 6 , the solidity J= c/s mayepe 
computed for various assumed values of b. By plotting these re- 
Sets, it is possible to determine b graphically for any choice 
of J . Table F I gives a set of such corresponding values for 
a p ee 45°. 

In the present example the solidity was arbitrarily chosen 
as unity. The parameters ei ae cascade geometry are now fixed. 

In order to obtain the point to point relation between the 
% and 53 plane it was necessary to reduce equation Fl to its 
real and imaginary parts. After a simple but somewhat lengthy 


calculation there is obtained the form, 


AS. 





“ (r+ %)cosQ@ —Zacosh7| +} {(% - r)sénQ} 


Fla 
= | oe | 
. \(r- V\cosQ +2Q sinh] +7 1% i) sin Of 
—+2tan® In ~ = 
(r-V%)cosO-Za sinh} +3 (% +r) sin@} 
tG e 
Where ea and pxC . 
THE STRAIGHT THROUGH FLOW 
The complex potential function necessary to produce a 
straight through flow with no circulation is 
a’ 
ur = D+2 eae In Sth + {n2t 7b F3 
er hs 5 =p $= Wp 


where vee represents the strength of the sources and sinks in 
the 5 plane. In the ¥ plane Vq is the velocity at infinity 
and s is the cascade spacing. The product VAS , therefore rep- 
sents the quantity of flow passing between two adjacent blades. 
The deviation of the velocity along any streamline from the 
value Va is a measure of the error in locating the bus bars a 
finite distance from the cascade. If the ratio ee along 
any streamline, at a position along the streamline corresponding 


to the location of the bus bars, is negligible then the error due 


to the bus bars' location is negligible. 


ae 





The stagnation streamline was chosen for investigation as 
it has the most radical changes in velocity and direction. Cal- 
culations for this streamline in the 5 plane occur where QG=O, 
fm and b >Yr2>Q . Setting Q= for convenience, and placa 
O=TT, the mapping equation reduces to: 


Pyescosie. | Ar + Yr) = 2cosht 
AT "| 2eoshY 4 (t+ A) : 


Fine 


uy = 5 sin 1, | {’- eS 
AT Z2sinht +r- YW) 


In order to facilitate visualization of perpendicular distances 
| 
from the cascade centerline a simple rotation of axes in the a 


plane is introduced, 


wautiy where x= X'cosB-ysinB 





F4 
and Y= x'sinB+y’cos6. 
4 
The complex velocity in the % plane is, 
ee: = 
Sur — = \y =< VE pane t ® F3a 
7 


mem the selections of O=TI, Q = 1.00 and 6 = 45° which deter- 


mines 8 for this case, 


Ve Pai + stanh T(t!) a 
VA < Om 4 tanh* 7)" 


The magnitude of this ratio does not change under a rotation 





of the axes. 
The value of Y in the § plane was varied from 1.0 to b. 


This gave corresponding points in the % or %y plane from 


Be2 





equations Flb and F4 respectively. At the same time the ra- 
tO aoe vas computed through the use of equation F3b. 

ra ee of */s and corresponding values of Na =Mry 
for the stagnation streamline of the cascade are — ia 
Table FII with Y as the argument. Fig. F3 provides a graphic 
display of the results. The cascade in the figure is drawn 
to scale. 

The working surface as finally evolved had bus bars located 
16 inches from the cascade centerline. Two sizes of profiles, 
6 and 12 inch chord, were used in the investigation. The 
solidity for the flat plate cascade has been selected as unity 


Maereby limiting the values of the ratio x/s to 16/12 and 16/6. 


These two cases are shown in Fig. F3. 


1 Ba, 





TABLE Fa 


CORRESPONDING VALUES OF b and 0 at p = 45° 


b ( at B = 49° 
2.00 6560 
1.95 ,6739 
1.90 6930 
1.85 ‘ailge 
1.80 7350 
1.75 7586 
1.70 eer 
1.65 8111 
1 5610 8413 
1.55 8737 
1.50 9100 
1.45 9506 
1.40 9966 
1.35 1.0491 
1.30 | lier 
1.25 1.1846 
1.20 1.2765 
1.15 1.3969 
1.10 1.5699 
1.05 1.8716 
1.00 00 
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TABEE Eas 


us @) 
CORRESPONDING VALUES OF = AND See at 6 = leoneme B = 45 


Via 
O02 nase oho, Dies 
1.04 tee OO! 8209 
OO +e OO e12un 
i. O8 +.3274 .6724 
iO #33594 0926 
m2 Ee ooL9 20324 
~— 3535 cascade edges 
1.14 a JOS8 .4793 
1 BR +.3794 4113 
8 +. 3948 BeOS 
20 +.4111 ReNUSS 
i 2 +.4345 72085 
24 +.4503 eee 
20 +.4732 . 1846 
IL APRS *.4994 2LOU5 
i oO OUD 1189 
eS 2 +.5684 Oy 
1.34 e201 90 0639 
50 arp eps) .0494 
38 oO i sOTGS 
it. 39 ora .0055 


ee 
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